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The  Escherichia  coli  RNA  degradosome  recognizes  and  degrades  RNA  through  the  coordination  of  four 
main  protein  components,  the  endonuclease  RNase  E,  the  exonuclease  PNPase,  the  RhIB  helicase  and 
the  metabolic  enzyme  enolase.  To  help  our  understanding  of  the  functions  of  the  RNA  degradosome, 
we  quantified  expression  changes  of  >2300  proteins  using  mass  spectrometry  based  shotgun 
proteomics  in  E.  coli  strains  deficient  in  rhIB,  eno,  pnp  (which  displays  temperature  sensitive  growth),  or 
rne(1-602)  which  encodes  a  C-terminal  truncation  mutant  of  RNase  E  and  is  deficient  in  degradosome 
assembly.  Global  protein  expression  changes  are  most  similar  between  the  pnp  and  rhIB  mutants, 
confirming  the  functional  relationship  between  the  genes.  We  observe  down-regulation  of  protein 
chaperones  including  GroEL  and  DnaK  (which  associate  with  the  degradosome),  a  decrease  in 
translation  related  proteins  in  A  pnp,  A  rhIB  and  rne{  1-602)  cells,  and  a  significant  increase  in  the 
abundance  of  aminoacyl-tRNA  synthetases.  Analysis  of  the  observed  proteomic  changes  points  to  a 
shared  motif,  CGCTGG,  that  may  be  associated  with  RNA  degradosome  targets.  Further,  our  data 
provide  information  on  the  expression  modulation  of  known  degradosome-associated  proteins,  such  as 
DeaD  and  RNase  G,  as  well  as  other  RNA  helicases  and  RNases  -  suggesting  or  confirming  functional 
complementarity  in  some  cases.  Taken  together,  our  results  emphasize  the  role  of  the  RNA 
degradosome  in  the  modulation  of  the  bacterial  proteome  and  provide  the  first  large-scale  proteomic 
description  of  the  response  to  perturbation  of  this  major  pathway  of  RNA  degradation. 


Introduction 

The  RNA  degradosome  is  a  vital  component  of  mRNA  degradation 
in  Escherichia  coli :1  and  a  main  contributor  to  post-transcriptional 
gene  regulation.  The  catalytic  core  enzyme  of  the  RNA  degradosome, 
the  endonuclease  RNase  E  (RNase  E),  is  conserved  amongst  proteo- 
bacteria,  but  homologs  have  also  been  identified  in  archaea  and 
plants.1  It  consists  of  an  N-terminal  catalytic  region  and  a  C-terminal 
non-catalytic  region  which  provides  a  scaffold  for  binding  of  other 
proteins,  RNA  or  other  cellular  components  (Fig.  1).  The  canonical 
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RNA  degradosome  in  E.  coli  is  a  500-700  kDa  complex  formed  by 
the  assembly  of  polynucleotide  phosphorylase  (PNPase,  the  product 
of  the  pnp  gene)  which  functions  as  an  exonuclease,  the  DEAD-box 
helicase  RhlB,  and  the  glycolytic  enzyme  enolase  (Eno)  on  the 
C-terminal  end  of  the  RNase  E.1  Cells  lacking  RNase  E  are  not 
viable,  most  likely  due  to  defects  in  tRNA  transcript  processing.2 
However,  truncation  mutants  lacking  the  C-terminal  scaffolding 
domain  and  thus  incapable  of  assembling  the  degradosome 
complex  are  viable  although  they  exhibit  slow  growth  and 
temperature  sensitivity.2  RNase  III  (Rnc)  and  RNase  G  (Rng) 
can  complement  rne  mutants  whereas  RNase  R  (Rnr)  and  RNase 
II  (Rnb)  can  substitute  for  PNPase.1 

Several  other  proteins  have  been  suggested  to  associate  with 
the  degradosome  or  play  a  role  in  RNA  degradation.  For 
example,  the  poly(A)  polymerase  PcnB  interacts  with  the 
DEAD-box  helicases  SrmB,  CsdA  and  RhlE.3  RNase  E  co-precipitates 
with  Ppk,  DnaK  and  GroEL  in  sub-stoichiometric  amounts  relative 
to  the  major  degradosome  components.1  Pull-down  experiments 
have  further  identified  a  number  of  proteins  bound  to  degrado¬ 
some  components  [e.g.  see  Fig.  3),  but  the  role  of  these  interactions 
is  unclear.4 
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Fig.  1  The  Escherichia  coli  RNA  degradosome  complex.  The  figure  shows  the  four  main  RNA  degradosome  components,  RNase  E,  Enolase,  RhIB,  and  PNPase.  RNase  E 
comprises  a  catalytic  N-terminal  region  and  a  non-catalytic  C-terminal  region  to  which  the  other  three  proteins,  RhIB,  Eno  and  Pnp,  bind.  In  addition  to  interacting  with 
RNase  E,  RhIB  and  Pnp  also  physically  interact  with  each  other.7  The  RNase  E  truncation  mutant  contains  only  the  catalytic  N-terminal  domain;  the  other  three  mutants 
are  gene  deletions. 


While  the  physiological  significance  of  the  E.  coli  degradosome 
is  well-established,  there  are  still  many  unresolved  questions 
regarding  its  function.1  For  example,  the  exact  role  of  the  RNA 
degradosome  in  modulating  mRNA  decay  and  thus  protein  abun¬ 
dance  in  vivo  is  unclear.5  The  significance  of  the  enolase  binding  to 
the  C-terminal  half  of  RNase  E  has  not  been  defined  apart  from  its 
role  in  modulating  the  mRNA  stability  of  the  glucose  transporter 
ptsG  during  metabolic  stress.6  Enolase  does  not  bind  RNA,  and 
only  5-10%  of  the  total  enolase  in  cell  extracts  is  complexed  with 
the  degradosome.5  RhlB  and  PNPase  can  also  form  a  complex  that 
degrades  double-stranded  RNA  independent  of  RNase  E.7  Further¬ 
more,  it  is  not  known  how  the  individual  component  enzymes 
cooperate  in  degrading  RNA  or  how  changes  in  the  degradosome 
composition  affect  the  decay  of  different  transcripts.8 

To  answer  some  of  these  questions  at  a  systems  level,  several 
studies  have  analyzed  changes  in  transcript  abundance  and 
mRNA  half-lives  in  E.  coli  strains  deficient  in  degradosome 
components,  namely  in  cells  carrying  null  mutations  in  pnp ,  eno  or 
rhlB  or  an  RNase  E  C-terminal  tmncation  that  abrogates  degrado¬ 
some  assembly  (rae(l-602)).9-12  Interference  with  degradosome 
function  may  lead  to  accumulation  of  un-degraded  RNAs  in  the 
cell  and  an  increase  in  their  RNA  half-lives.  However,  transcript 
studies  observed  a  wide  range  in  responses,  with  a  substantial 
number  of  transcripts  that  decreased  in  half-lives  or  did  not 
change  in  concentration,11  suggesting  that  the  regulation  of 
RNA  degradation  is  complex. 

The  ultimate  effect  of  mutations  that  influence  mRNA  half- 
life  is  a  change  in  the  steady  state  proteome.  Here  we  present  a 


large-scale  quantitative  analysis  of  changes  in  the  proteome  of 
strains  lacking  the  canonical  degradosome  components  and 
compare  our  data  with  the  published  information  on  the  mRNA 
half-lives  in  these  mutants.  Deep  protein  abundance  analysis 
was  carried  out  using  high-resolution  mass  spectrometry  and  a 
label-free  quantitation  method  developed  in  our  lab13,14  to 
estimate  the  concentration  changes  of  >2300  proteins.  These 
data  provide  novel  information  on  how  the  function  of  the 
degradosome  shapes  the  E.  coli  proteome. 

Methods 

Data  collection 

E.  coli  strains  carrying  pnp ,  rhlB ,  eno  mutations  or  rne{ 1-602) 
were  kindly  provided  by  Dr  S.  N.  Cohen  (Stanford  University)10 
(Table  1).  Cells  were  grown  at  30  °C  in  M9  media  supplemented 
with  0.2%  tryptone,  0.2%  glycerol,  1  mM  MgS04  and  0.0001% 
(w/v)  thiamine.  40  mM  succinate  was  added  to  strains  K10  and 
DF261  as  described  in  ref.  10.  Cell  pellets  were  harvested  when 
the  culture  absorbance  reached  A600  =  0.6,  the  cells  were 
resuspended  in  lysis  buffer  (25  mM  Tris  HCl,  pH  7.5,  1  mM 
EDTA,  2.5  mM  DTT,  and  lx  Roche  (Nutley,  NJ)  Protease 
inhibitor)  and  then  lysed  by  passing  three  times  through  a 
French  Press  at  20  000  psi.  Soluble  fractions  were  collected  by 
centrifugation  (10  000  xg,  10  min,  at  4  °C)  and  diluted  to  4  mg  ml-1. 
The  protein  concentration  was  determined  by  Nanodrop  (Thermo 
Scientific).  Protein  samples  were  heat  denatured  at  95  °C  for  15  min. 
After  cooling  down  to  room  temperature,  trypsin  (Sigma)  was  added 


Table  1  Strains  carrying  mutations  in  the  major  degradosome  components.  All  experiments  were  conducted  as  described  in  ref.  10.  The  protein  expression  levels 
were  quantified  in  mutant  and  wild-type  strains  separately,  and  the  Z-score  calculation  from  ref.  14  used  to  estimate  significance  and  direction  of  the  expression 
changes.  See  ESI  for  assessment  of  data  quality 


Name  of 
experiment 

Mutant 

strain 

Genotype 

No.  of  proteins  quantified  in 
expression  change  (5%  FDR) 

Wild-type 

Genotype 

Aeno 

DF261 

K10  except  for  eno-2 

1576 

K10 

garBlO,  fhuA22,  ompF627(T2R),  fadL70l(T2R), 
relAl,  pit-10,  spoTl,  rrnB-2,  mcrBl,  creC510 

Apnp 

YHC012 

N3433  except  for  Tn5::pnp 

1651 

N3433 

lacZ,  relA,  spotl,  thil 

ArhlB 

SU02 

N3433  except  for  A  rhlB 

1540 

N3433 

lacZ,  relA,  spotl,  thil 

Truncated  rne 

BZ453 

SH3208  except  for  RNase  E 
truncation 

1545 

SH3208 

his  AtrpE5(7) 
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at  1  :  50  (w/w).  Digestions  were  carried  out  at  37  °C  for  24  hours. 
The  sample  was  again  lyophilized  to  10  pi,  resuspended  in 
120  pi  buffer  (95%  H20,  5%  acetonitrile,  0.1%  formic  acid)  and 
filtered  through  a  Microcon-10  filter  at  12  000g.  The  samples 
were  stored  at  -80  °C  until  LC-MS/MS  analysis. 

Samples  (10  pi)  were  injected  into  an  LTQ-Orbitrap  Classic 
(Thermo  Electron)  mass  spectrometer  and  peptides  were 
resolved  using  a  5  to  90%  acetonitrile  gradient  over  five  hours 
via  reverse  phase  chromatography  on  a  Zorbax  BioBasic-18 
column  150  mm  x  0.10  mm  ID  (Thermo  Fisher).  Each  of  the 
runs  was  analyzed  independently  with  Bioworks  (Thermo 
Fisher),  searching  a  database  of  E.  coli  protein  sequences 
imported  from  the  NCBI  genomic  database.  The  results  were 
combined  for  analysis  by  PeptideProphet,15  ProteinProphet16 
and  post-processed  in  the  APEX  pipeline13,14  to  estimate  absolute 
and  differential  protein  expression  based  on  adjusted  spectral 
counts.  We  accepted  proteins  as  confidently  identified  if  the 
ProteinProphet  probability  was  above  a  cutoff  corresponding  to 
<5%  global  false  discovery  rate.  Relative  protein  expression 
changes  were  calculated  with  respect  to  measurements  in  the 
respective  wild-type  strain  (Table  1),  using  a  Z-score  calculation 
described  by  Lu  et  al.14  The  Z-score  adjusts  the  error  model  to 
the  absolute  concentration  (spectral  counts).  A  Z-score  of  |Z|  > 
1.96  corresponds  to  a  P-value  <  0.05.  Note  that  Z-scores  can  be 
calculated  even  if  one  of  the  samples  (mutant  or  wild-type)  has 
zero  counts  measured,  but  a  fold-change  in  expression  cannot. 
We  collected  at  least  five  technical  replicates  per  sample,  i.e. 
repeat  mass  spectrometry  measurements  which  were  pooled  for 
the  quantitative  analysis.  These  replicates  also  included  repeat 
sample  preparations.  Replicate  concentration  measurements 
correlated  well  (R2  =  0.74  to  0.94,  not  shown)  and  were  pooled  for 
further  analysis.  The  raw  data  are  available  at  Tranche  (Proteome- 
Commons),  data  hash:  zO+qZFSpP8Z5XzXrSSTm3KYpF2fE2J- 
8F3pDUZpSL5WtSOdIdQ+lF/eGXJHUoXMsRtcdkHzfFmStufDON- 
4S9BaOSjQqOAAAAAAAADIA==.  More  information  on  parameters 
of  the  mass  spectometry  experiment  and  a  discussion  of  data 
quality  is  provided  in  Section  1  (ESI|). 


Computational  analysis 

The  mRNA  microarray  expression  data  and  estimated  half-lives 
were  obtained  from  Bernstein  et  al.10  and  consist  of  relative 
mRNA  expression  values  for  1825  genes  with  missing  values 
estimated  using  a  K  nearest  neighbor  method.  Tiling-array 
data  for  rne{ 1-602)  cells  were  from  Stead  et  al.11  and  consist 
of  2720  genes  for  which  significant  up-  or  down-regulation  is 
reported. 

All  data  were  analyzed  with  a  combination  of  Perl  scripts,  R,  and 
the  software  Perseus.17  The  combined  Z-score  for  a  protein’s  expres¬ 
sion  change  (Fig.  3  and  4)  was  calculated  as  the  sum  of  all  Z-scores 
across  all  four  mutants  divided  by  the  square  root  of  the  number  of 
Z-scores  available  for  the  particular  protein:  Zcomb  =  [Z  °f all 
contributing  mutants]  /  [number  of  contributing  mutants] .  The  com¬ 
bined  Z-score  is  a  measure  for  the  general  trend  in  the  data:  we 
assume  that  a  protein  (mRNA)  which  is  consistently  up-regulated 
across  the  four  RNA  degradosome  mutants  is  more  likely  a  direct 
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Fig.  2  Protein  and  RNA  expression  response  of  degradosome  mutants.  The 
heatmap  shows  the  Z-scores  of  the  protein  expression  data  (first  four  columns)  of 
the  four  degradosome  mutants  mapped  to  tiling-array  based  mRNA  expression 
data  of  the  RNase  E  truncation  mutant1 1  (rightmost  column).  Significant  function 
enrichments  (adjusted  P-value  <  0.01,  compared  to  entire  E.  coli  genome)  are 
indicated  for  the  four  largest  clusters  (yellow  boxes),  with  cluster  sizes  in  brackets. 
Only  genes  without  missing  data  are  shown  ( N  =  904).  Colors  saturate  at  values 
of  |2|  to  illustrate  statistically  significant  expression  changes  (P-value  <  0.05)  - 
actual  Z-scores  range  from  -15  to  8.  Fold-changes  of  expression  range  from  -6.0 
to  5.7  on  a  log  scale  (base  2).  t-rne  denotes  the  truncation  mutant  me(1-602). 


target  of  the  degradosome  than  proteins  (mRNAs)  with  low  com¬ 
bined  Z-scores. 

For  clustering  of  the  data  in  Fig.  2,  up  to  two  missing  values 
were  accepted  for  the  proteomics  data,  and  all  missing  values 
in  the  tiling  array  data  were  substituted  by  0.  Data  were 
re-normalized  to  a  mean  of  0  and  standard  deviation  1.  Clustering 
was  based  on  the  Euclidean  distance  measure  using  an  average 
linkage  algorithm  (Perseus17),  with  a  1.8  Euclidian  distance  cutoff 
to  determine  genes  clusters  with  similar  expression  profiles.  Gene 
function  analysis  was  performed  for  the  largest  clusters  using 
FuncAssociate.18  Sequence  motifs  were  analyzed  using  MEME19 
and  fuzznuc.20 


Results 

Effects  of  A pnp,  Aeno ,  A rhlB  or  rne( 1-602)  on  the  E.  coli 
proteome  -  data  quality  and  general  results 

We  grew  E.  coli  strains  carrying  Aeno ,  A  rhlB,  Apnp ,  or  rne{ 1-602) 
alleles  to  the  mid-exponential  phase  and  used  mass  spectro¬ 
metry-based  proteomics  of  cell  lysates  to  quantify  soluble 
protein  expression  changes  relative  to  the  respective  parental 
(wild-type)  strains,  in  a  manner  similar  to  the  published  mRNA 
analysis.12  We  provide  statistical  scores  of  expression  changes 
for  a  total  of  2390  proteins,  quantifying  >1500  proteins  per 
mutant  (Table  1).  Fig.  2  shows  the  direction  and  statistical 
significance  (Z-scores)  of  the  protein  expression  changes 


752  |  Mol.  BioSyst.,  2013,  9,750-757 


This  journal  is  ©The  Royal  Society  of  Chemistry  2013 


Published  on  24  January  2013.  Downloaded  by  University  of  Texas  Libraries  on  28/08/2013  20:44:44. 


View  Article  Online 


Paper 


Molecular  BioSystems 


gN 

e  & 

Q.  c 


sf5 

fl 

5>  a 

CO  o 


200 

1.00 

0.00 

-I.OQ^ 


Aeno 


.  i.. 


a  r 


s; 

s 

fl 

C 

g 

<75 


1.00 

000 


•'V7 


..n.,  1,1  -  JJ 

r  M  ^  ^  ^  ^  f  &  &  «H  cf 

i  1 1 :  -At  I  ■ 


Apnp 


■m 
1.00 
000 
■'0%f 

m 

» ti 

1 1  100 

*  000 

200 


2  ~ 

S  & 

§sS 

if  * 

l| 

c/>  5 


Sf  J 

32 

If 


Arhlb  |" 

ll 

.  1 

1 

<?  * 

' 

'  ^  1 

1  <r  </ 

j 

H 

|  <S*  ^  ^ 

truncated 


>oo  _ 

.  [  ■ 

>nn  ® 


;ated  me  I 

ll  .1  I  ll 

^  s  </  M  #  ^  ^ 


A  A 


Fig.  3  Expression  changes  of  RNA  degradosome  proteins  and  related  helicases  and  RNases.  The  figure  shows  the  significance  (Z-score)  and  direction  of  the  protein 
expression  change  of  the  four  main  mutants  as  determined  by  mass  spectrometry  experiments  for  the  degradosome  components  and  other  known  RNases  and  RNA 
helicases  that  can  potentially  buffer  for  the  loss  of  gene  function.  Color  codes  follow  the  designation  in  Fig.  1 .  A  Z-score  of  |Z|  >  1 .96  is  significant  within  a  P-value  of 
0.05.  The  figure  shows  a  range  of  Z  between  {-2,  2}  for  clarity  reasons;  actual  Z-scores  range  from  - 1 5  to  8  (ESI t).  Corresponding  fold-changes  of  protein  expression 
can  be  found  in  Fig.  S4  (ESIt).  The  value  for  the  gene  in  the  respective  mutant  is  left  out  (see  Fig.  SI,  ESI  t).  Triangles  indicate  the  significance  of  the  combined  Z-score  at 
the  5%  level  (see  Methods). 


between  the  mutant  strains  and  their  respective  wild  type  cells 
for  the  subset  of  904  proteins  that  are  completely  characterized 
(no  missing  values,  capped  at  \Z\  =2  corresponding  to  a 
P-value  =  0.05). 

Label-free  quantitative  proteomics  reports  protein  expression 
changes  with  a  precision  of  about  2-  to  3-fold.14  The  proteomics 
experiment  conducted  here  correctly  detects  expression  changes 
of  the  four  mutant  proteins  with  highly  significant  Z-scores:  as 
expected,  the  proteins  encoded  by  eno,  pnp,  and  rhlB  are  absent 
in  the  respective  null  strains  whereas  in  rne{ 1-602)  cells  the 
truncated  RNase  E(l-602)  protein  is  over-expressed,  consistent 
with  previous  observations12  (Fig.  SI,  ESI|).  For  further  quality 
control,  we  conducted  western  blot  experiments  on  two  soluble 
proteins,  GlmS  and  CspE  (Fig.  S2,  ESI|).  For  both  proteins,  the 
proteomics  and  western  blot  data  are  in  strong  agreement  for  the 
PNPase  mutant,  and  less  so  for  the  RhlB  mutant.  GlmS  concen¬ 
tration  drastically  increases  in  the  PNPase  mutant  compared  to 
wild-type  both  in  the  proteomics  and  western  blot  data.  This 
increase  contrasts  the  concentration  decrease  at  the  mRNA  level 
which  can  be  explained  by  strong  regulation  of  GlmS  translation 


and  mRNA  degradation  by  two  small  non-coding  RNAs  (GlmY 
and  GlmZ)  which  in  turn  interact  with  the  RNA  degradosome.21 
This  result  is  also  consistent  with  our  recent  study  in  which  we 
searched  for  mutants  which  can  up-regulate  rraB  transcription. 
In  our  study,  we  isolated  a  glmS  mutant  in  which  a  Tn5 
insertion  disrupts  the  coding  region  of  the  GlmS  protein.22 
RraB  disrupts  PNPase  associated  with  RNA  degradosome  both 
in  vivo23  and  in  vitro  (L.  Zhou,  G.  Georgiou,  unpublished).  The 
proteomics  data  are  consistent  with  the  notion  that  there  is 
negative  feedback  between  PNPase  and  GlmS. 

In  general,  the  fold-change  in  expression  values  between 
mutant  and  wild-type  for  a  protein  appears  to  be  of  similar, 
if  not  larger,  magnitude  to  that  of  the  corresponding  mRNA 
(Fig.  S4  and  S6,  EST\).  Examples  of  fold-changes  of  protein 
expression  for  some  proteins  are  provided  in  Fig.  S5  (ESIfl. 
About  25%  of  the  proteins  detected  an  increase  in  relative 
abundance  in  all  four  mutant  strains  (A eno,  ArhlB ,  A pnp, 
rne{ 1-602))  (Fig.  2).  When  grouping  protein  expression  profiles 
according  to  the  Euclidean  distance,  the  A  rhlB  and  A pnp  cells 
proved  to  be  more  similar  to  each  other  than  any  other  mutant 
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Fig.  4  Expression  changes  of  proteins  associated  with  the  RNA  degradosome.  The  figure  shows  the  significance  (Z-score)  and  direction  of  the  protein  expression 
change  of  several  genes  with  known  association  with  the  RNA  degradosome  or  RNA  degradation.  A  Z-score  of  |Z|  >  1.96  is  significant  within  a  P-value  of  0.05.  The 
figure  shows  a  range  of  Z  between  {-2,  2};  actual  Z-scores  range  from  -1 5  to  8  (ESIt).  Triangles  indicate  the  significance  of  the  combined  Z-score  at  the  5%  level  (see 
Methods). 


pairs  (Fig.  2;  Fig.  S4,  ESIf),  supporting  the  known  physical  and 
functional  interactions  between  the  two  enzymes.7 

Functional  analysis  of  the  four  large  clusters  in  Fig.  2 
indicates  a  global  down-regulation  of  translation-related  proteins 
(chaperones,  ribosomal  subunits),  primarily  in  the  Apnp  and  A rhlB 
mutants  (cluster  B).  Such  down-regulation  is  consistent  with 
regulatory  coupling  between  mRNA  availability  and  translation. 
However,  the  proteomics  data  also  demonstrate  an  up-regulation  of 
14  of  the  20  E.  coli  aminoacyl-tRNA  synthetases  in  all  four  mutants 
(cluster  A,  Fig.  2;  Fig.  S3,  ESI,|  P-value  <  0.001).  To  the  best  of 
our  knowledge,  a  possible  link  between  the  RNA  degradosome 
and  aminoacylation  of  tRNAs  has  not  yet  been  described.  We 
hypothesize  that  this  up-regulation  may  occur  as  the  cell’s 
effort  to  counteract  the  imbalanced  concentrations  of  mRNAs 
and  mature  tRNAs  due  to  impaired  degradosome  function,  or 
possibly  as  part  of  a  general  stress  response  that  produces 
alarmones.24 

Next,  we  examined  likely  degradosome  targets  for  enrichment 
of  sequence  motifs  which  may  be  recognized  by  components  of 
the  RNA  degradation  machinery.  We  derived  these  likely  targets 
by  combining  the  Z-scores  describing  the  expression  change  in 


the  four  individual  mutants  into  one  combined  Z-score.  Across 
the  protein  and  RNA  expression  data,  as  well  as  mRNA  half-life 
data,25  we  consistently  observe  enrichment  in  the  CGCTGG 
motif  in  which  the  center  CTG  triplet  is  the  most  conserved  part 
(datafile  3,  ESIt).  This  motif  occurs  more  frequently  in  mRNAs  of 
short  half-life  than  mRNAs  with  long  half-life,  and  it  is  enriched 
in  sequences  that  increase  in  either  mRNA  or  protein  expression 
levels  when  RNA  degradosome  function  is  impaired.  Thus,  the 
CGCTGG  can  be  loosely  associated  with  direct  targets  of  the  RNA 
degradosome.  The  ESI|  discusses  this  observation  and  analysis 
of  other  sequence  features. 

The  expression  response  in  degradosome  components  and 
protein  family  members 

Fig.  3  indicates  the  protein  expression  changes  for  the  four 
main  degradosome  components  as  well  as  for  other  E.  coli 
RNases  and  RNA  helicases  detected  by  the  proteomics  analysis. 
Note  that  not  all  RNases  and  RNA  helicases  are  associated  with 
the  degradosome,  and  that  the  figure  caps  expression  changes 
at  a  Z-score  of  >2  which  corresponds  to  significant  expression 
changes  at  a  P-value  of  <0.05. 
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RNase  R  (Rnr)  is  a  3'  to  5'  exoribonuclease  closely  related  to 
RNase  II  (Rnb)  which  is  important  for  decay  of  mRNA  with 
extensive  secondary  structures.26  Since  the  coordination 
between  RNA  degradosome  components  organized  by  the 
C-terminal  region  of  RNase  E  is  required  for  decay  of  structured 
RNA  substrates,  it  is  tempting  to  assume  that  the  observed 
increase  in  the  cellular  concentration  of  RNase  R  in  cells 
expressing  the  truncated  RNase  E(l-602)  (Fig.  3)  is  due  to  a 
defect  in  structured  RNA  observed  in  this  strain.  Similarly, 
RNase  BN  (Rbn)  is  implicated  in  the  processing  of  tRNA 
precursors.27  RNase  BN  abundance  decreases  in  both  the  rhlB 
and  pnp  mutants  (Fig.  3),  but  is  unchanged  in  rne{ 1-602)  (or  in 
the  eno  mutant).  The  RNase  BN’s  (Rbn)  inverse  expression 
pattern  with  RNase  II  (Rnb)  is  consistent  with  their  overlapping 
functions  in  vivo  28  down-regulation  of  one  enzyme  may  be  com¬ 
pensated  for  by  up-regulation  of  the  other  enzyme  and  vice  versa. 

E.  coli  has  five  known  DeaD-box  RNA  helicases  (ribonu- 
cleases):  RhlB,  RhlE,  SmrB,  DeaD  (also  called  CspA),  and 
DbpA.1  Null  mutants  of  these  helicases  result  in  viable  cells, 
but  can  affect  the  growth  rate.29  DeaD,  SrmB  and  RhlE  have 
been  reported  to  associate  with  RNase  E  in  vitro.5,30,31  DeaD  and 
RhlE  have  evidence  for  their  role  in  RNA  degradation;1  they  can 
complement  the  RNA  helicase  function  in  the  absence  of 
RhlB,32  whereas  SrmB  cannot.33  The  binding  sites  for  DeaD, 
SrmB  and  RhlE  within  the  C-terminal  end  of  RNase  E  are 
different  from  those  of  RhlB,  and  overlap  with  the  enolase 
binding  site  (residues  833-850,  Fig.  1).  We  observe  an  overall 
increase  in  the  protein  expression  levels  of  DeaD,  SrmB,  and 
RhlE  across  the  four  mutants  which  is  largely  consistent  with 
the  functional  complementarity  between  the  enzymes  and 
degradosome  components  (Fig.  3). 

HrpA  is  also  an  RNA  helicase  (although  not  from  the  DEAD- 
box  family)  that  functions  in  mRNA  processing.34  Its  expression 
level  is  drastically  down-regulated  in  both  the  rhlB  and  pnp 
mutants,  similar  to  the  effect  observed  for  another  RNA  heli¬ 
case,  YejH  (Fig.  3).  The  expression  patterns  of  these  two 
helicases  contrast  those  of  the  DeaD-box  helicases. 

RNase  III  (Rnc)  is  involved  in  ribosomal  RNA  processing,35  an 
activity  which  overlaps  with  the  essential  role  of  RNase  E  in  E.  coli. 
Similar  to  RNase  R,  we  observe  significant  up-regulation  of  RNase 
III  in  the  rae(l-602)  mutant  consistent  with  their  complementary 
roles  in  rRNA  processing  and  in  ribosome  biogenesis  (Fig.  3). 
Although  RNase  G  (Rng)  is  an  RNase  E  parolog  that  cannot  fully 
complement  RNase  E  function,36  RNase  G  concentration  also 
increases  in  three  of  the  four  mutants  (Fig.  3).  Both  the  endonu¬ 
clease  RNase  G  and  the  DeaD  helicase  are  up-regulated  in  the  eno, 
pnp  and  me  mutants,  a  finding  consistent  with  a  recent  genetic 
analysis  by  Tamura  et  al.  showing  that  mutation  in  DeaD  can 
suppress  RNase  E  in  a  temperature-sensitive  me  mutant.12 

RNase  A  ( rna )  is  involved  in  tRNA  and  rRNA  processing,37 
and  our  data  show  slight  down-regulation  in  the  RNase 
E  mutant  (Fig.  3).  The  oligonucleotide-ribonuclease  encoded 
by  orn  acts  downstream  of  the  RNA  degradosome  where  it 
regenerates  single  nucleotides  from  oligonucleotides.38  Its 
expression  levels  changes  are  similar  to  RNase  E  and  RNase  G 
except  for  the  rne  mutant  (Fig.  3). 
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The  expression  response  of  other  degradosome  associated 
proteins 

We  also  assembled  a  number  of  genes  which  have  some  known 
association  with  the  RNA  degradosome  based  on  pull-down 
interaction  data  and  other  evidence,1  and  for  which  our  proteomics 
data  provided  information  (Fig.  4).  Since  these  proteins  are  neither 
RNases  nor  RNA  helicases,  which  are  discussed  above,  their  role  in 
RNA  degradation  is  presumably  indirect. 

In  particular,  three  chaperones  (GroEL,  GroES,  and  DnaK) 
and  two  ribosomal  proteins  (RplB,  RpsC)  are  strongly  down- 
regulated  in  the  four  mutants  (somewhat  less  in  the  RNase  E 
truncation)  (Fig.  4)  -  corresponding  to  cluster  B  in  Fig.  2.  The 
RNA  polymerase  RpoA  is  significantly  up-regulated  in  the  rhlB 
and  pnp  mutants,  but  down-regulated  in  the  eno  and  rne{ 1-602) 
mutant  strains  (P-value  <  0.05)  (Fig.  4),  suggesting  feedback 
between  transcription  and  degradation  regulation.  The  poly(A) 
polymerase  PcnB  is  up-regulated  in  all  four  mutants  in  our  data 
as  might  be  expected  since  these  strains  display  varying  degrees 
of  defects  in  RNA  processing  and  PcnB  may  counteract 
this  defect.  The  addition  of  poly(A)  tails  (by  PcnB)  generally 
destabilizes  RNA  in  bacteria,39  and  it  facilitates  degradation  of 
REP-stabilizers,  i.e.  structural  elements  preventing  degradation, 
with  the  help  of  RhlB.1 

The  protein  expression  response  to  perturbation  of  the  RNA 
degradosome  differs  from  the  RNA  response 

We  compared  the  protein  expression  changes  to  published 
steady  state  transcriptome  and  tiling  array  data10,11  to  investi¬ 
gate  the  relationship  between  transcript  and  protein  expres¬ 
sion.  Since  partially  degraded  RNA  can  still  hybridize  to  arrays, 
and  translation  efficiency  of  partially  degraded  RNAs  is  likely 
affected  in  the  mutants  as  well,  we  did  not  expect  much 
agreement  between  the  protein  and  mRNA  data.  Indeed,  the 
protein  expression  data  do  not  correlate  with  either  mRNA 
expression  changes  or  mRNA  half-life  changes  measured  by 
microarrays  (Fig.  S4,  ESIfl.10 

mRNAs  that  are  not  processed  in  the  rne{ 1-602)  mutant 
(which  is  unable  to  form  the  degradosome  complex)  may  still 
be  translated  into  proteins.  Inspection  of  the  published  tiling 
array  data11  in  Fig.  2  (rightmost  column)  indicates  that  about 
one  third  of  the  transcripts  for  which  protein  abundance  data 
could  be  obtained  are  possibly  processed  in  a  degradosome- 
dependent  manner  (Fig.  2,  clusters  A,  D,  (227  +  98)/904  =  36%): 
their  expression  levels  increase  in  tiling  array  in  the  rne{ 1-602) 
mutant.  Of  these,  a  large  number  of  transcripts  also  show 
positive  protein  expression  changes  in  the  rne{ 1-602)  mutant 
(Fig.  2,  second  column  from  the  right),  but  the  agreement 
between  tiling  array  data  and  protein  expression  data  for  the 
RNase  E  mutant  is  generally  low.  Only  one  quarter  (227/904, 
cluster  A)  of  the  proteins  in  the  dataset  are  largely  up-regulated 
in  their  protein  expression  patterns  across  all  four  mutants  and 
in  the  tiling  array  data.  The  remaining  set  of  proteins  may 
accumulate  at  the  mRNA  level,  but  their  functional  half-life 
does  not  seem  to  be  affected,  and  protein  expression  levels  are 
not  up-regulated  (cluster  D).  Similarly,  a  number  of  genes  are 


This  journal  is  ©  The  Royal  Society  of  Chemistry  2013 


Mol.  BioSyst.,  2013,  9,750-757  |  755 


Published  on  24  January  2013.  Downloaded  by  University  of  Texas  Libraries  on  28/08/2013  20:44:44. 


View  Article  Online 


Molecular  BioSystems 

either  unaffected  in  their  RNA  expression  levels  by  the  RNase  E 
truncation  or  down-regulated  (clusters  B  and  C),  and  these 
genes  may  be  degraded  by  pathways  other  than  the  canonical 
degradosome.  Again,  the  protein  data  match  the  mRNA  data 
only  partially  (cluster  B). 


Conclusions 

We  present  a  large-scale  proteomics  analysis  of  mutants  of  the 
E.  coli  RNA  degradosome,  complementing  existing  data  on 
mRNA  expression  and  half-life  changes.10  Since  measurements 
of  mRNA  half-lives  do  not  inform  on  translation  activity  of  the 
transcripts  and  the  final  effect  of  RNA  degradation  on  protein 
abundance,  proteomics  analysis  provides  an  important  additional 
angle  to  our  understanding  of  RNA  degradation.  The  proteome-wide 
data  provided  here  indicate  several  trends  and  novel  associations 
that  can  be  followed  up  by  future  experiments.  Expression  data 
of  >2300  proteins  reveal  substantial  differences  between  the 
effects  of  perturbed  RNA  degradation  on  transcript  and  on 
protein  concentrations:  fewer  than  half  of  the  genes  in  our 
dataset  appear  to  follow  the  canonical  model  of  RNA  degrada¬ 
tion  (Fig.  2).  A  trivial  explanation  of  this  discrepancy  lies  in 
technical  measurement  noise.  However,  the  estimated  error  in 
protein  expression  measurements  is  smaller  than  the  observed 
variation  between  experiment  and  control  here,14,40  supporting 
the  existence  of  true  biological  trends. 

The  complexity  of  E.  coli  RNA  degradation  can  be  partially 
explained  by  the  up-regulation  of  several  RNA  helicases  and 
RNases  that  may  buffer  for  loss  of  enzyme  function  in  the 
degradosome  mutants.  We  are  able  to  provide  proteomic 
characterization  for  RNA  helicases  and  RNases  known  to  have 
such  buffering  roles  (DeaD,  Rnr,  Rng),1  but  also  suggest  possible 
complementary  functionality  for  additional  enzymes  (HrpA,  YejH, 
Rna,  Orn)  -  but  their  precise  interaction  with  the  degradosome 
needs  to  be  tested  in  future  experiments.  We  provide  two  novel 
findings  that  have  not  yet  been  discussed  in  the  context  of 
the  RNA-degradosome.  First,  we  observe  a  large  number  of 
aminoacyl-tRNA  synthetases  to  be  up-regulated  in  the  pnp ,  rhlB , 
and  rne  mutants,  suggesting  a  link  to  tRNA  metabolism.  Second, 
we  observe  enrichment  of  the  CGCTGG  motif  in  sequences  that 
are  likely  direct  targets  of  the  RNA  degradosome  -  the  role  of 
which  is  entirely  unknown. 

Our  results  provide  insights  into  the  cellular  role  and 
structure  of  the  RNA  degradosome  itself.  Since  only  a  small 
fraction  of  the  PNPase  and  enolase  molecules  participate  in  the 
degradosome  complex,5,41  it  is  unsurprising  that  their  mutant 
expression  profiles  are  highly  discordant  with  those  of  the  other 
mutants,  and  highlight  the  roles  of  these  proteins  outside  the 
degradosome.  However,  the  profiles  of  the  rhlB  and  pnp 
mutants  correlate  better  with  each  other  than  the  profiles  of 
the  other  mutants  (R  =  0.72),  consistent  with  their  physical  and 
functional  interaction  independent  of  the  RNase  E  scaffold.7 
Interestingly,  this  correlation  between  the  rhlB  and  pnp  mutant  is 
observed  for  protein  expression  (this  study)  and  mRNA  half-lives,10 
but  not  for  published  mRNA  expression  data  (Fig.  S4,  ESI|), 
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suggesting  a  highly  dynamic  and  complex  system  of  RNA 
degradation. 

The  lack  of  correlation  amongst  the  datasets  also  suggests 
extensive  regulation  at  the  level  of  translation  and  protein 
degradation  as  well  as  RNA  degradosome  functions  outside 
the  removal  of  mRNAs.  Indeed,  even  in  unperturbed  wild-type 
bacterial  systems,  a  low  correlation  between  mRNA  and  protein 
concentrations  has  indicated  that  translation  and  protein 
degradation  regulation  account  for  substantial  amounts  of  protein 
expression  variation.14,42  The  same  is  true  for  eukaryotes; 
transcription  and  RNA  degradation  account  for  only  half  of 
the  variation  in  steady-state  protein  expression  levels.40,43 
Further,  mRNA  half-lives  and  transcript  concentrations  are 
inversely  correlated  in  E.  coli,9  suggesting  additional  layers  of 
regulatory  complexity. 

Finally,  the  discordance  between  protein  and  RNA  expression 
levels  observed  is  also  consistent  with  a  model  that  distinguishes 
between  the  chemical  and  functional  half-life  of  mRNAs.1  The 
mRNA  expression  data  indicate  changes  in  the  chemical  half-life, 
i.e.  the  physical  presence  of  the  transcript.  The  functional  half- 
life  may  be  very  different  from  the  chemical  half-life,  affecting 
the  translation  activity  associated  with  the  mRNA  and  influen¬ 
cing  the  final  protein  expression  levels.  An  mRNA  may  be 
physically  present  in  the  cell,  but  translationally  inactive  and 
thus  non-functional.  Our  proteomics  data  deliver  supporting 
evidence  for  the  distinction  between  chemical  and  functional 
half-lives  of  mRNAs,  but  detailed  biochemical  studies  will  be 
required  in  order  to  confirm  and  further  characterize  this  mode 
of  regulation. 
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